Long-range correlated random field and random anisotropy O(N) models: A functional 

renormalization group study 
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We study the long-distance behavior of the O(N) model in the presence of random fields and random 
anisotropies correlated as ~ l/x d ~" for large separation x using the functional renormalization group. We 
compute the fixed points and analyze their regions of stability within a double e = d — 4 and a expansion. We 
find that the long-range disorder correlator remains analytic but generates short-range disorder whose correlator 
develops the usual cusp. This allows us to obtain the phase diagrams in (d, a, N) parameter space and compute 
the critical exponents to first order in e and a. We show that the standard renormalization group methods with a 
finite number of couplings used in previous studies of systems with long-range correlated random fields fail to 
capture all critical properties. We argue that our results may be relevant to the behavior of 3 He-A in aerogel. 

PACS numbers: 64.60.Ak, 64.60.Fr, 75.10.Nr, 74.25.Qt 



I. INTRODUCTION 



The effect of weak quenched disorder on large-scale prop- 
erties and phase diagrams of many-body systems attracted 
considerable attention for decades. Despite significant efforts, 
there still remain many open questions. The prominent exam- 
ple is the O(N) model in which an TV-component order pa- 
rameter (the magnetization in the spin system notation) is cou- 
pled to a symmetry-breaking random field. For N = 1, it is 
known as the random field Ising model (RFIM). 1 For N > 1, 
one has to distinguish the random field (RF) case where the 
order parameter couples linearly to disorder, and the random 
anisotropy (RA) case, where the coupling to disorder is bilin- 
ear. These models are relevant for a variety of physical sys- 
tems such as amorphous magnets, 2 diluted antiferromagnets 
in a uniform external magnetic field,? liquid crystals in porous 
media, 4 critical fluids in aerogels, 5 nematic elastomers^ and 
vortex phases of impure superconductors. 7 For N = 1, the RA 
model reduces to the random-temperature model, where the 
randomness couples to the local energy density as, for exam- 
ple, in diluted ferromagnets. 8 In contrast to the systems with 
random-temperature-like disorder, the RF and RA models suf- 
fer from the so-called dimensional reduction (DR). A straight- 
forward analysis of the Feynman diagrams giving the leading 
singularities yields to all orders that the critical behavior of 
the RF O(N) model in d dimension is the same as that of the 
pure system in d — 2 dimensions. 9 Consequently, the lower 
critical dimension is (N = 1) = 3 for Ising-like systems 
and df c R (N > 1) = 4 for systems with continuous symmetry. 
This can elegantly be demonstrated using supersymmetry. 10 
However, simple Imry-Ma arguments show that the lower crit- 
ical dimension of the RFIM is d lc (N = 1) = 2^ The de- 
viation from DR is also confirmed by the high-temperature 
expansion. 12 Thus DR breaks down, rendering standard field 
theoretic methods useless. 

Another known problem where the perturbation theory is 
spoiled by DR is elastic manifolds in disordered media. There, 
two methods were developed to overcome difficulties re- 
lated to DR: the Gaussian variational approximation (GVA) in 



replica space and the functional renormalization group (FRG). 
The GVA is supposed to be exact in the limit N — ► oo. 
Unfortunately, this approach when applied to the RF prob- 
lem leads to very complicated equations which do not al- 
low us to compute the critical exponents. 13 Considering the 
RF O(N) model, Fisher 14 showed that expansion around the 
lower critical dimension d\ c = 4 generates an infinite num- 
ber of relevant operators which can be parametrized by a 
single function. However, he found that the corresponding 
one-loop FRG equation has no analytic fixed point (FP) solu- 
tion. Only recently, using the progress in the elastic manifold 
problem j 15 ' 16,17 ' 18 it was realized that the scaling properties 
of systems exhibiting metastability are encoded in a nonan- 
alytic FR Feldman 19 has shown that, indeed, in d — 4 + e 
and for N > N c « 3 (a more precise computation 20 gives 
N c = 2.834 74), there is a nonanalytic FP with a cusp at the 
origin. This FP provides the description for the ferromagnetic- 
paramagnetic phase transition in the RF O(N) model and 
allows one to compute the critical exponents which are dif- 
ferent from the DR prediction. Recently, Le Doussal and 
Wiese 20 extended the FRG analysis to two-loop order. The 
extension beyond one-loop order is highly nontrivial due to 
the nonanalytic character of the renormalized effective action, 
which leads to anomalous terms in the FRG equation^ The 
two-loop calculations were also independently performed in 
Ref.HT], and the truncated exact FRG was proposed in Ref.l22i 
The more accurate analysis of the FRG flows for the RF model 
showed that for N > N* = 18 + 0(e) there is a crossover 
to a FP with weaker nonanalyticity resulting in the DR critical 



exponents 
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A similar picture was found for the RA 



0(N) model with the main difference that N c = 9.4412 and 

N* = ooM 

A more peculiar issue concerns the phase diagram of the 
RF and RA models below d\ c . It is known that for the 
RF model and models with isotropic distributions of random 
anisotropies true long-range order is forbidden below d\ c = 4 
(for anisotropic distributions, long-range order can occur even 
below di c )^ However, quasi-long-range order (QLRO) with 
zero order parameter and an infinite correlation length can per- 
sist even for d < d\ c . For instance, the GVA predicts that 
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the vortex lattice in disordered type-II superconductors can 
form the so-called Bragg glass exhibiting slow logarithmic 
growth of displacements. 25 This system can be mapped onto 
the RF 0(2) model, in which the Bragg glass corresponds to 
the QLRO phase. Indeed, for N < N c and d < d\ c , the 
FRG equations have attractive FPs which describe the QLRO 
phases of RF and RA models. 26 In order to study the tran- 
sition between the QLRO phase and the disordered phase ex- 
pected in the limit of strong disorder, one has to go beyond the 
one-loop approximation. The truncated exact FRG (Ref. 22) 
and the two-loop FRG (Ref. l20l) performed using a double 
expansion in •v/fij and N — N c give access to a different 
singly unstable FP which is expected to control the transi- 
tion. Both methods give qualitatively similar pictures of the 
FRG flows: the critical and attractive FPs merge in some di- 
mension dj* c (N) < d\ c that is therefore the lower critical di- 
mension of the QLRO-disordered transition. For the RF 0(2) 
model, both methods predict approximately the same lower 
critical dimension d* c 3.8(1), and thus, suggest that there is 
no Bragg glass phase in d = 3. However, one has to take cau- 
tion when extrapolating results obtained for small \/\s\ and 
N — N c . Moreover, in contrast to the model of Refs. 2q and 
Eol which belongs to the so-called "hard-spin" models, the sys- 
tem studied in Ref. l22l corresponds to "soft spins", and thus, is 
expected to belong to a different universality class. In terms of 
vortices, the soft spin RF model allows for topological defects 
which destroy the Bragg glass. 

Most studies of the RF and RA models are restricted to 
either short-range (SR) correlated disorder or uncorrected 
pointlike defects. However, real systems often contain long- 
range (LR) correlated disorder or extended defects in the 
form of linear dislocations, planar grain boundaries, three- 
dimensional cavities, etc. Systems with anisotropic orienta- 
tion of extended defects can be described by a model in which 
all defects are strongly correlated in Ed dimensions and ran- 
domly distributed over the remaining d — Ed dimensions. The 
case Ed = is associated with uncorrelated pointlike defects, 
while extended columnar or planar defects are related to the 
cases Ed = 1 and 2, respectively. The critical behavior of the 
O(N) model with random-temperature-like extended defects 
was studied in Refs. l27U28 1 29 30 using a perturbative RG anal- 
ysis in conjunction with a double expansion in e = 4 — d and 

In the case of an isotropic distribution of disorder, power- 
law correlations are the simplest example with the possibil- 
ity for a scaling behavior with new FPs and new critical 
exponents. The critical behavior of systems with random- 
temperature disorder correlated as \/x d ~ a for large separa- 
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tion x was studied in Refs. 31 32.33. The power law can be as- 
cribed to extended defects of internal dimension a, distributed 
with random orientation in d dimensional space. In general, 
one would probably not expect a pure power-law decay of cor- 
relations. However, if the correlations of defects arise from 
different sources with a broad distribution of characteristic 
length scales, one can expect that the resulting correlations 
will, over several decades, be approximated by an effective 
power law. 31 Power-law correlations with a noninteger value 
a can be found in systems containing defects with fractal di- 



mension a. For example, the behavior of He in aerogels is 
argued to be described by an XY model with LR correlated 
defects. 35 This is closely related to the behavior of nematic 
liquid crystals enclosed in a single pore of aerosil gel which 
was recently studied in Ref. |36|, using the approximation in 
which the pore hull is considered a disconnected fractal. The 
FRG was used to investigate the statics and dynamics of elas- 
tic manifolds in media with LR correlated disorder in Ref.l37l 
The critical behavior of the 0(N) model with LR correlated 
RF was studied in Refs. [38ll39ll40[ However, the methods used 
in these works fail to describe properly the case of SR cor- 
related RF. Therefore, there is a necessity to reexamine the 
critical behavior of the LR RF model using methods which 
are successful in the SR case. 

In the present paper, we study the LR correlated RF and RA 
O(N) models using the FRG to one-loop order. The paper 
is organized as follows. In Sec. [TTJ we introduce the LR RF 
and RA models, and derive the FRG equations. In Sec. [TTTJ 
we study the LR RF model. In Sec. IIV1 we consider the LR 
RA model and discuss the application to superfluid 3 He in 
aerogels. The final section summarizes our results. 



II. MODEL AND FRG EQUATIONS 

The large-scale behavior of the O(N) symmetric spin sys- 
tems at low temperatures can be described by the nonlinear 
model with the Hamiltonian 



H[S]= J d d x 



V(x, s) 



(1) 



where s(x) is the A^-component classical spin with a fixed- 
length constraint s 2 = 1. V(x, s) is the random disorder po- 
tential, which can be expanded in spin variables as follows: 



V{x,s) 



E -^\{x) Sll (x)...s H Xx). (2) 



The corresponding coefficients have simple physical inter- 



pretation: is a random field, h\^' is a random second- 
rank anisotropy, and h^' are general /ith tensor anisotropies. 
As was shown in Ref. [lj, even if the system has only fi- 
nite number of nonzero the RG transformations gen- 
erate an infinite set of high-rank anisotropies preserving the 
symmetry with respect to rotation s — > — s if it is present 
in the bare model. For instance, starting with only second- 
rank anisotropy corresponding to the RA model, all even-rank 
anisotropies will be generated by the RG flow. We will reserve 
the notation RA for the systems which have this symmetry and 
the notation RF for the systems which do not. In the present 
work, we consider the case of Gaussian distributed long-range 
correlated disorder with zero mean and cumulants given by 



(2) 



h\tiM<U( x ') = ^5 lin ...5^y^5{x~x>) 



» 



00; 



(3) 



with g (x — x 1 ) ~ l/|x — x'\ d <T . For the sake of convenience, 
we fix the constant in Fourier space, taking g(q) = l/q° . 
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To average over disorder, we introduce n replicas of the 
original system and compute the replicated Hamiltonian 



Tin 



d d x 



~ E( w ° ) 2 ~ 2^ E MM*) ■ M*)) 

a a,b 



J d d x g{x- x')R 2 (s a {x) ■ s b {x')) 

a.b 



,(4) 



where Ri(z) — YluSi Z>J "- ^he P ro P ert i es °f the original 
disordered system (fl]) and (0 can be extracted in the limit 
n — > 0. According to the above definition of the RF and RA 
models, the functions Ri(z) are arbitrary in the case of the RF 
and even for the RA. Power counting suggests that d\ c = 4+ a 
is the lower critical dimension for both models4£ 

At criticality or in the QLRO phase, the correlation func- 
tions of the order parameter exhibit scaling behavior. In con- 
trast to the models with temperature-like disorder in the mod- 
els under consideration, the connected and disconnected cor- 
relation functions may scale with different exponents. This re- 
flects metastability and the breaking of the DR. For instance, 
the connected two-point function behaves as 



(%) • s(-q)) (%)) • (s(-q)) 



while the disconnected function scales as 



(%)) • (s(-q)) (%)) • (s(-q)) 



7 -2+r) 



-4+J7 



(5) 



(6) 



Here, s(q) is the Fourier component of the order parame- 
ter and the angle brackets stand for the thermal averaging. 
Schwartz and Soffer 41 proved that the exponents of the RF 
model obey the inequality 2r/ > fj. Since in the RA case the 
coupling to disorder is bilinear, the Schwartz-Soffer inequal- 
ity cannot be applied directly to r\ and fj. In the RA case, it is 
convenient to introduce the correlation functions of the form 
Eqs. (O and $6$ not for s, but for the field = [si] 2 — 1 /N, 
and define exponents 772 and 772 which satisfy the Schwartz- 
Soffer type inequality: 2r/2 > f]2- 2b Vojta and Schreiber— gen- 
eralized this inequality to correlated RF and obtained a more 
restrictive bound 2-q — a > fj since a > 0. For RA models 
similar arguments lead to 2r/2 — a > 772- 

To derive the one-loop FRG equations, we straightfor- 
wardly generalize the methods developed in Refs.QjlOand 
l37l to model (0). We express the order parameter s a as a 
combination^ 6 . 



s a (x) = n a (x) y/l- 7r a 2 (z) + 7? a (x) 



(7) 



of a fast field if a fluctuating at small scales A < q < Aq which 
is orthogonal to a slow field n a of unit length, changing at 
scales q < A. Here, Ao is the UV cutoff and A <C Ao. The 
field fi a can be considered as the coarse-grained order param- 
eter (local magnetization) whose fluctuations at low tempera- 
ture are small, (tt 2 ) <C 1. Integrating out the fast variables 7? a , 
we rescale in such a way that the effective Hamiltonian of the 
slow fields ri a would have the structure of the bare Hamilto- 
nian It is convenient to change variable to z = cos 0. The 



FRG equations to first order in e and u are given by 

dtRitf) = -eRM) + \[R'{^) + i? 2 '(0)] 2 - AR'{{<j>) 



(8a) 



- {N - 2) 1 2 ARi {(f>) + AR[ (0) cot <p 
1 -[R'M+R'MY 



2 sir/ 



diR 2 (4>) - -(e-a)R 2 (<t>)-{(N-2)[2R 2 (cj)) 

+R 2 (4>)cot(t>] +R 2 \(p)\A, (8b) 

where di :— —d/d In A. We have absorbed the factor of 
l/(87r 2 ) in redefinition of R and introduced 



A = i?'/(0) +i? 2 '(0). 



(9) 



In terms of the variable <fi, the functions Ri (<fi) become peri- 
odic with period 2ir in the RF case and tt in the RA case. The 
flow equation for the temperature to one-loop order reads 



dihxT 



-(d-2) - (N — 2) A. 



(10) 



According to Eq. (ITOD . the temperature is irrelevant for d > 2 
and sufficiently small A. Although we expect A — 0(e, S) in 
the vicinity of a FP, one has to take caution whether the found 
FP survives in three dimensions. 26 The scaling behavior of 
the system is controlled by a zero-temperature FP of Eqs. dSat 
and d8b]l [R*,R£,A*], such that d e R* = 0. An attractive FP 
describes a phase, while a singly (unidirectionally) unstable 
FP describes the critical behavior. The critical exponents are 
determined by the FRG flow in the vicinity of the FP and to 
one-loop order are given by 

r] = -A*, f} = -E- (N- 1)A*, (11a) 

772 = -(N + 2)A* , 772 = -£-2A^*. (lib) 

It is convenient to introduce the following reduced variables: 

r,($ = Ri{4>)/{e-a), (12a) 
a = A/(s-a), (12b) 
e = e/(e-a). (12c) 

To check the stability of the FP [r*,r|,a], we linearize the 
flow equations around this FP: ri((f>) = r*(<f>) + yi((f>) and 
obtain 

X yi (4>) = -e Vl {<t>) + [rficf) + r* 2 "(<f>)] [v'M + vZW] 
-ay'i{<t>) - a rl"{4>) — (N — 2){2a r^) 
+a r*' (cj)) cot <p + 2ayi (</>) + ay[ ((/>) cot 4> 
-[rJ'W+rS'W] 

x[yi(0)+y 2 (0)]/sin 2 0}, (13a) 
A2/2W = -y 2 {c}>)-a Q {{N-2)[2rl{ct>) + r%{<f>)cot<t>} 
+r* 2 "(<t>)} - a{(N- 2) [2y 2 (cf>) +M) cot <£\ 
+2/ 2 '(0)}, (13b) 
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TABLE I: LR RF model above the lower critical dimension d\ c — 
4 + a. The FP values of ?"*"(0), r*"(7r) and the relevant eigenvalue 
Ai computed numerically for different N and £. The last column is 
the relevant eigenvalue A^ obtained from the truncated RG scheme 
of Ref. 39 and computed using Eq. l !16b . 



N 


£ 


rf'(O) 




Ai 


A? 


4 


1.271 


-1.0000 


0.5811 


1.271 


2.429 




2 


-0.4657 


0.0891 


1.218 


2.000 




3 


-0.3320 


0.0041 


1.198 


1.618 




4 


-0.2668 


-0.0207 


1.192 


1.414 


5 


2 


-0.1743 


0.0192 


1.167 


1.366 




3 


-0.1132 


-0.0073 


1.160 


1.225 




4 


-0.0819 


-0.0138 


1.144 


1.158 


6 


2 


-0.0941 


0.0080 


1.145 


1.215 




3 


-0.0549 


-0.0058 


1.127 


1.135 



1.30 



1.25 



1.20 



£ 115 



1.10 



1.05 



1.00 



a 
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12 14 

N 



16 



20 22 



where we have introduced the eigenvalues A (measured in 
units of £ - a) and denned a = y'{(0) + y£(0). The FP 
is attractive if all A,; fulfill the inequality (e — a)\ < 0. 
A singly unstable FP has only one eigenvalue Ai such that 
(e — a)Xi > 0, which determines the third independent expo- 
nent 



u = l/[Xi(e-a)]. 



(14) 



This exponent characterizes the divergence of the correlation 
length in the vicinity of transition. 



III. LONG-RANGE RANDOM FIELD O(N) MODEL 

We now focus on the phase diagram and critical behavior 
of the LR RF 0(N) model. Kardar et al^ studied the critical 
behavior of the LR RF 0(N) model using a e = d uc — d 
expansion around the upper critical dimension d uc = 6 + <r. 
They found that, to lowest order in e, the critical properties 
are that of a pure system in d — 2 — a dimensions; however, 
this generalized DR was found to fail at higher orders. Chang 
and Abrahams 39 applied to the LR RF O(N) model a low- 
temperature version of the RG. Expanding around the lower 
critical dimension d\ c = 4+<j, they obtained the RG recursion 
relations for the three parameters T, Ai, and A2, which in 
our notation correspond to Aj = — R"(0). As we have shown 
in the previous section, the truncated RG neglects an infinite 
number of relevant operators. They found the nontrivial zero- 
temperature FP, which in our notation Si = A.; /(e — a) reads 



61 



1 



(e-l)(N-3y 



i(N -3) - N + 2 
(e-l)(N-3)* ■ 



(15) 



The correlation length exponent v is determined by Eq. ( fT"4"i >, 
with the relevant eigenvalue Ai given by 



N- 2 
N-3 




-2-i(N- 3)] 
e 2 (AT-3) 2 ' 



(16) 



Note that the expression for the relevant eigenvalue reported 
in Ref. 39 is incorrect and gives values which are several times 



FIG. 1 : The stability regions of various FPs corresponding to differ- 
ent patterns of the critical behavior above the lower critical dimen- 
sion, e > a. The borderline between the SR and LR FPs is com- 
puted numerically using Eq. ( I22t ; the borderline between the SR TT 
and LR TT regions is given by Eq. J23l > and that between the LR and 
LR TT regions by Eq. l !25t . Inset: Schematic phase diagram on the 
(e, cr) plane for a particular value of iV € (3, 18). The solid line is 
a — e; the dashed and dotted lines are given by Eqs. | |22| | and l !25t . 
respectively. 



larger than that computed using Eq. ( fl~6*l ). Although the trun- 
cated RG scheme when applied to the model with only SR cor- 
related RF (A 2 = 0) results in the DR FP A^ DR = e/(N— 2) 
and exponent ^dr = 1/e, the exponent ( fl~6l ) differs from the 
generalized DR prediction t'DR = l/(e — <r). Nevertheless, 
one may doubt about the applicability of the truncated RG to 
the LR RF O(N) model even if the DR is broken. 

We now reexamine the long-distance behavior of the LR 
RF O(N) model by means of the full one-loop FRG derived 
in Sec.[ll] Equation fl8bl is linear in the function r2(<f>) and 
can be solved analytically. The FP solution fulfilling the RF 
boundary condition, i.e., 2-7T periodicity, is given by 



r5(0=-r2"(O)cos& 



(17) 



which is an analytic function. Note that the analyticity of the 
LR part of the disorder correlator was also revealed in the elas- 
tic manifold problem. 37 This also gives us the FP value of a, 



1 



*LRRF 



AT-3' 



(18) 



which, following Eqs. ( II 1 al l, completely fixes the values of 
the critical exponents 77 and fj, 



e — a 
N-3'' 



2e- (N — l)o 
N-3 ' 



(19) 



Exponents ( fT9l ) satisfy the generalized Schwartz-Soffer in- 
equality at equality. This is at variance with the SR RF mod- 
els, where the Schwartz-Soffer inequality was found to be 
strict, 19 but in agreement with the results for the LR RF spheri- 
cal model with long-range interactions. 42 One may conjecture 
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that the generalized Schwartz-Soffer inequality is satisfied as 
equality for any N including the LR RF Ising model^ but 
there is no argument that this persists in higher orders of loop 
expansion. 

Let us first discuss the critical behavior of the LR RF model 
above the lower critical dimension, e > a. Note that — Olrrf 
exactly coincides with 6 1 + 5% given by Eqs. (Tl~5b - Therefore 
the truncated RG would give the same values of 77 and fj as the 
full FRG, at least to first order in e and a, if these exponents 
would be computed in Ref. |39[ To obtain the function r*(0) 
and the amplitude ^"(O), we integrate Eq. ( |8ab numerically. 
Since the coefficients of Eq. (I5at are singular at 4> = and ir, 
to compute the solution in the vicinity of these points, we use 
expansions of r\(4>) in powers of \<j>\ and (tt — <fi) 2 , respec- 
tively. The expansion around is completely determined by 
the value of r*"(0), 



± 



2(N 2 - 47V + 2>)rl" (0) + N - 1 <'(O)0 2 
2(N - 3)[{N - 3)i - 2(N - 2)] + 2 
y/rt"(0)(£-l){N- 3) 2 + l 



-H 3 + O(0 4 ),(2O) 

while to get the explicit expansion around tt we need know 

r*"(0) and r*"(ir): 

rt{4>) = {N-l) 

{ H"(0) + (JV- 3) + I}' + 2r 1 *"M(AT - 3) 



rr'W(7T 



2(iV - 3)[{N - 3)e - 2(2V - 2)] 

^+O[(7T-0) 4 ]. 



(21) 

From Eq. d20l i. we see that the SR part of the disorder cor- 
relator is nonanalytic at small <f>, so that we have to distin- 
guish the left and right derivatives. In what follows, we adopt 
r[ n) (0) = r[ n) (0+). We use numerical integration to continue 
the solutions given by expansions ( f20b and OH inside the in- 
terval [0, tt] and match them by adjusting the shooting param- 
eters r*"(0) and r*"(ir). Only the series with "+" in Eq. (f2Qb 
can be matched with the solution computed using expansion 
( 1211 . Following the third term of Eq. d20l l, the FP solution 
exists only if — r*"(0) < 5*. We found that this inequality is 
always strict, however, the difference r*"(0) — (— 5*) becomes 
smaller in the limit of large N. Thus the truncated RG fails to 
give the correct values of r*"(0) and ^"(0) although it gives 
the correct sum. Consequently, the relevant eigenvalue Ai and 
the exponent v obtained from the truncated RG are also ex- 
pected to be incorrect. The computed values of r*"(0) and 
r*"(7r) are shown in TableU] 

Let us check the stability of FPs. First, we examine the 
stability of the SR FP [r^ R ,r 2 = 0,a^ R = r^ R (0)] found 
numerically in Ref. |26|. Linearized about the FP, Eq. d 1 3bl > 
can be solved analytically giving 3/2 = cos and A = — 1 — 
o,g R (iV — 3). The SR FP is stable against the introduction 
of LR disorder if (e — er)A < 0. Taking into account that 
VSR = -(e-cr)ag R andfysR = -e-(N-l)(e-cr)a? 



SR' 



can rewrite the criterion of the SR FP stability in the following 
form 



Disordered 




FIG. 2: Phase diagram of the RF model below the lower critical di- 
mension, d < 4 + a (s < a). 



which was derived using general RG scaling arguments in 
Ref. H^. The regions of stability of the LR and SR FPs are 
shown in Fig. Q] Tissier and Tarjus (TT) argued that for 
N > 18, the cuspy SR FP becomes more than once unstable, 
but a different critical (singly unstable) SR TT FP arises with 
the function (4>) being only p ^ N times differentiable at 
the origin. 21 Although this weak nonanalyticity should reflect 
some metastability in the system, 20 it leads to the DR critical 
exponents. According to Eq. ( f22b . the SR TT FP is stable with 
respect to the LR correlated disorder for N > 18 and 



Ke<(JV-2)/(JV-3). 



(23) 



We now check the stability of the LR RF FP 
[rj (0), r^O), a = -1/(N - 3)]. Substituting the LR 
FP in Eqs. (1 1 3ab and ( I13bb . we obtain y2{4>) — cos (j> and 



A 



a = 



l + r*"(0)(iV-3)' 



(24) 



To compute the eigenfunction yi((j>) and eigenvalue A, we 
solve Eq. ( 1 1 3at numerically using shooting and imposing the 
27r-periodic boundary condition. The obtained values of Ai 
are shown in Table |U As one can see from the table, the rele- 
vant eigenvalues Ai computed using the full FRG, and conse- 
quently, the exponent v differ significantly from that obtained 
using the truncated RG. The critical exponents are expected 
to be continuous functions of e, a, and N. In the region 
controlled by the SR FP the correlation length exponent is 
^sr = l/ e ' independent on NM- The first line in TableU cor- 
responds to a point on the borderline separating the regions of 
stability (see Fig. Q]), and thus, on the borderline we indeed 
have vsr = ^lr- We checked that the second eigenvalue A2 
computed at the LR FP vanishes exactly on the borderline, in- 
dicating that the LR FP becomes twice unstable in the region 
dominated by the SR FP. We also looked for an analog of the 
TT phenomena in the presence of LR correlated disorder and 
found that for 



a < 2?ysR - i7sr, 



(22) 



e > 



2VN + 7 + > 
N-3 



(25) 
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there is a crossover to a new singly unstable LR FP of a TT 
type with the function R\ ((f)) being only p ~ N times differ- 
entiable at the origin and R\ {<t>) ~ cos 4>- The values of R'{(0) 
and i? 2 '(0) at the LR TT FP are simply given by Eq. < TT3T >. 
Thus, the LR TT FP leads to the critical exponents ( [T9| > and 
the relevant eigenvalue (TToT i obtained from the truncated RG. 

Finally, we discuss the LR RF O(N) model below the lower 
critical dimension, e < a. For N < 3, the FRG equations 
have two attractive FPs: the SR and LR, which describe the 
SR QLRO and LR QLRO phases, respectively. The phase di- 
agram computed from the stability analysis of different FPs is 
depicted in Fig. |2] The exponents characterizing the power- 
law behavior of the connected and disconnected correlation 
functions in the LR QLRO phase are given by Eq. ( fT9l ). The 
physically interesting case N = 2 describing the Bragg glass 
phase in the presence of the LR correlated disorder was con- 
sidered by one of the authors in Ref. H3- In the Bragg glass, 
the displacements u(x) of a periodic structure (e.g., vortex 
lattice) grow logarithmically as (u(x) — u(0)) 2 = .Ad In |x|. 
Note that in Ref. [37| the period was fixed to 1 while in the 
present work it is 2tt, so that the relation between the univer- 
sal amplitude Ad and the exponent 77 is given by 77 = 2ir 2 Ad- 
For N = 2, the crossover from SR QLRO with the critical 
exponents t]sr = \e\ir 2 /9 and fjsn = |e|(l + tt 2 /9) to LR 
QLRO with the exponents i%r = \e\ + a and ?7lr = 2|e| + a 
happens for e < 9/ir 2 (see Fig. [2j». 



IV. LONG-RANGE RANDOM ANISOTROPY 0(N) 
MODEL 



In this section, we study the LR RA case which corresponds 
to 7r-periodic functions ri((f>). Equation dFbl has a family of 
(at most) tt -periodic solutions which can be expressed as poly- 
nomials in cos <j) of power p with 



1 



Ap 2 + 2(p-l)(N -2) : 



P=l,2, 



For instance, forp = 1, we have 



1 



1], 



1 

4' 



(26) 

(27) 
(28) 



a LRRA 

and forp = 2, we have — 1/(2N + 12) and 

(2) 2(iV + 6)rp ),/ (0)-l 

r2 " 8(N +1)(N + 2)(N + 6) 

x [3 - 6(N + 2) cos 2 <f> + (N + 2)(N + 4) cos 4 </>]. 

The stability analysis shows that due to the inequality a* < 
a* = Olrra f° r P ^ 2 and N > 2, all FPs with p > 2 are 
unstable. It can be easily seen for N — 2 when the functions 
T^-f 1 (p > 2) become (7r/p)-periodic, and thus, are unstable 
with respect to a 7r-periodic perturbation. Henceforth we con- 
sider only the LR FP determined by Eqs. ( f27T > and (f28j), which 
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-0.5 
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LR QLRO 
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FIG. 3: Phase diagram of the RA model below the lower critical 
dimension, d < 4 + a (e < er). 



give the following values of the critical exponents 



?7LR 



?72LR 



—7— . ^LR = 

(a-e)(N + 2) 



-(N- 1) - -(JV + 3), (29) 



N 



V2LR = -e 



-(tr-e).(30) 



Exponents ( f30b satisfy the generalized Schwartz-Soffer in- 
equality 2ri2 — cr > fj-2 also at equality. Analogously to 
the LR RF model, in order to obtain the function r\{(j)), we 
solve Eq. dial ) numerically using shooting and imposing the u- 
periodic boundary condition. Since the coefficients of Eq. ( fSab 
are singular at <f> = 0, we use an expansion of r\ ((f)) in powers 
of 1 0| which reads 



(N — 1)[1 — 8r!"(0)] rl"{0)<j) 2 



16(JV - 2 + 2e) 



±A 1 + 



16r*"(0)(e- 1) |# 



AT • 



12 



O(0 4 ). (31) 



As one can see from Eq. (l3TT >, the SR disorder correlator r I ((f)) 
has a cusp at the origin. Only the solution with "+" in Eq. (|3T1 l 
fulfills all conditions. 

We now check the stability of the SR RA FP with respect 
to the LR correlated disorder. Analogously to the RF case, 
linearized around the FP, Eq. (II 3bb allows for an analytical 
solution giving yi = N cos 2 0—1 and A = — 1 + 4agR- The 
SR FP is stable if (e — <r)X < 0, which can be rewritten as 



(32) 



CT < 27/2SR - »/2SR- 



Above the lower critical dimension, e > a, inequality 
holds for all N > N c — 9.4412 so that the SR FP is stable 
with respect to the weak LR correlated disorder. Although the 
LR correlated disorder does not change the critical behavior 
for N > N c , it can modify the critical behavior which should 
exist for N < N c but which is not accessible in the one-loop 
approximation. 1 "^^ 

Below d\ c , the LR RA model cannot develop a true long- 
range order, but there can exist different types of QLRO. 
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The SR QLRO phase is controlled by the SR RA FR with 
r*(<f>) = rg R ((f>) and r% (4>) = computed in Ref. |26j. It was 
found to be stable in the subspace of the SR correlated disor- 
der for N < N c — 9.4412. We find that it is also stable with 
respect to the LR correlated disorder for a fulfilling inequality 
02t . The new LR QLRO phase is controlled by an attractive 
LR RA FP with 7r-periodic functions r\ ^ and 7-| ^ 0, 
which we computed numerically by integrating Eq. d8ab with 
the initial condition given by Eq. ( l3TT > and using r*" (0) as a 
shooting parameter. The exponents describing the correlation 
functions in the LR QLRO phase are given by Eqs. d29l and 
(f30b . Note that below d\ c we have rj > 0. To check the stabil- 
ity of the LR RA FR we substitute it in Eq. d 1 3bb and obtain 
2/2 = N cos 2 4> — 1 and 



2XN 



«o 



4r*"(0) - 1' 



(33) 



The eigenfunction yi(4>) and eigenvalue A are computed nu- 
merically using shooting. This allows us to determine the sta- 
bility regions for different QLRO phases, which are shown in 
Fig. [3] As one can see from Fig. [3] in contrast to SR QLRO, 
LR QLRO can exist even for 4 < d < 4 + a and N > N c . 

The above results may be relevant for the behavior of 3 He in 
aerogels. It has recently been observed in NMR experiments 
that the phase A of 3 He confined in aerogel exhibits two dif- 
ferent types of magnetic behavior called c and / states4i De- 
pending on the cooling, one can obtain either a nearly pure c 
state or a mixed f + c state, which gives two overlapping lines 
c and / in the transverse NMR spectrum. Although the pure / 
state has not been observed, there is an evidence that the f + c 
state is inhomogeneous and consists of regions with two dif- 
ferent magnetic orders c and J. 44 The order parameter of the 
3 He-A can be parametrized by the complex vector if) and the 
real unit vector d, which characterize the orbital and magnetic 
anisotropy, respectively. 45 Only the orbital part of the order 
parameter given by the real unit vector 1 = (i/2)[if),xp*] inter- 
acts with the aerogel matrix, which can be treated as quenched 
random anisotropy disorder. The spin part d is not coupled 
directly to the disorder, but there is a weak spin-orbit (dipole) 
interaction between 1 and d which generates the NMR fre- 
quency shift. In Ref. 46, the existence of the two states was 
interpreted in terms of different "random textures" of the field 
1 in the Larkin-Imry-Ma state (QLRO phase in our notation). 
The measured dependencies of the transverse NMR signal on 
the tipping angle can be explained if one assumes that in the 
c state the vectors 1 and d are almost uncorrelated, while in 
the / state they are partially locked. This is possible if in the 
c state the characteristic length of the texture, i.e., the Larkin 
length, is much smaller than the characteristic length of the 
dipole interaction, L <C while in the / state they are 
of the same order. 46 However, the nature of different random 
textures exhibiting different characteristic length scales is not 
clear. Alternatively, one can try to interpret the / state as a 
network of topological defects pinned by the aerog el.i 6 - We 
argue that these two different states may be the SR and LR 
QLROs found above in the LR RA model. Indeed, aerogel 
is a porous medium formed by tangled silicon strands which 



exhibit a fractal mass distribution. NMR and small-angle X- 
ray scattering (SAXS) experiments give the mass-to-distance 
relation M ~ x df , with the fractal dimensions df varying 
in the range of 1.4 - 2.4 depending on the fabrication pro- 
cess. This scaling holds up to the fractal correlation length, 
which can exceed 100 nm. 47 Thus, the effective RA disorder 
is expected to be long-range correlated with a = dt up to the 
scale of the fractal correlation length or even more^ As a re- 
sult, the SR QLRO phase is unstable to formation of islands 
with LR QLRO (see Fig. [3]). To compute the Larkin length, 
one has to solve the flow equations ( T8at and d8bl starting from 
a particular bare disorder correlator and looking for the scale 
at which the cusp is developing. 37 We can estimate the Larkin 
length using Flory-type arguments. 26 Assuming that J is an 
elastic constant and R is a strength of disorder, we obtain 
L S k ~ (J/R) 2/E and L LR - (J/R) 2 ^ 6 -^ for the SR and 
LR disorders, respectively. For aerogel, we have e = 1 and 
a « 1.4 ... 2.4 so that the Larkin lengths in both phases may 
differ significantly. This can explain the experimentally ob- 
served coexistence of regions with different spin states. Fur- 
ther investigations, however, are clearly needed. 



V. SUMMARY 

In this work, we investigated the long distance properties of 
the O(N) model with random fields and random anisotropies 
correlated as l/x d ~ a for large separation x. We derived the 
functional renormalization group equations to one-loop order, 
which allow us to describe the scaling behavior of the models 
below and above the lower critical dimension d\ c = 4 + a. 
Using a double e = d — 4 and a expansion, we obtained the 
phase diagrams and computed the critical exponents to first 
order in e and a. For the LR RF model, we found that the 
truncated RG developed in Ref. 39 to study the critical behav- 
ior above the lower critical dimension is able to give the cor- 
rect one-loop values of exponents ?/ and fj, but not the phase 
diagram and the critical exponent v except for the region con- 
trolled by the weakly nonanalytic LR TT FR Thus, although 
the truncated RG overcomes the dimensional reduction, it fails 
to reproduce all properties which can be obtained using the 
functional renormalization group. We found a new LR QLRO 
phase existing in the LR RF model below the lower critical 
dimension for N < 3 and determined the regions of its sta- 
bility in the (e, er, N) parameter space. We obtained that the 
weak LR correlated disorder does not change the critical be- 
havior of the RA model above d\ c for N > N c — 9.4412, but 
can create a new LR QLRO phase below d\ c . The existence 
of two QLRO phases in LR RA systems may explain the two 
different states of 3 He-A in aerogel observed recently in NMR 
experiments. 44 However, many questions are still open. In 
particular, the paramagnetic-ferromagnetic transition should 
exist also for N < N c , though it was not found in the one- 
loop approximation. Even for the SR correlated disorder, it is 
still unclear whether it remains perturb ative42i2£ It would also 
be interesting to find a connection with the replica symmetry 
breaking pictured 
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